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Abstract 

Wc further develop the numerical algorithm for computing the gauge connection 
of slope-stable holomorphic vector bundles on Calabi-Yau manifolds. In particu- 
lar, recent work on the generalized Donaldson algorithm is extended to bundles 
with Kahler cone substriicturc on manifolds with /i^'^ > 1. Since the computation 
depends only on a one-dimensional ray in the Kahler moduli space, it can probe 
slope-stability regardless of the size of h^'^. Suitably normalized error measures 
are introduced to quantitatively compare results for different directions in Kahler 
moduli space. A significantly improved numerical integration procedure based on 
adaptive refinements is described and implemented. Finally, an efficient numerical 
check is proposed for determining whether or not a vector bundle is slope-stable 
without computing its full connection. 
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1 Introduction 

In this paper, we explore Af = 1 supersymmetric vacua of Es x Eg heterotic string [Il[2] 
and M-theory |3l [U [5] [6] [T] [8] . The four-dimensional effective theory is specified by a 
Calabi-Yau threefold X and a slope-stable holomorphic vector bundle, V. The detailed 
structure of the low energy theory is determined [9] by the choice of a Ricci-flat metric, 
g, on the threefold and an = 1 supersymmetry gauge connection. A, on the vector 
bundle. Existence proofs, such as Yau's theorem [lUj for the Ricci-flat metric and the 
Donaldson-Uhlenbeck-Yau theorem |1H [T2] for the Hermitian Yang-Mills connection, 
provide us with numerous examples of such geometries. However, the explicit metric 
and gauge connection are not known analytically, except in very special cases |13l [T^ 
[151 [m [T71|18] . The difficulty of determining these quantities has presented an obstacle 
to the systematic search for realistic heterotic vacua. Even for a known vacuum, this 
has precluded the computation of physically relevant parameters in the effective theory, 
such as the Yukawa couplings. 

In recent years, the development of sophisticated numerical approximation schemes 
have provided a new approach to these problems [T9l[20t[2Tl[22l[251[2ll[25t[26l[27l[28]. 
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With the development of powerful new algorithms and modern computer speed, it is 
now possible to numerically approximate Ricci-flat metrics and Hermitian Yang-Mills 
connections to a high degree of accuracy. We will refer to these tools collectively as the 
"generalized Donaldson algorithm" . Using them, the structure of the four-dimensional 
effective theory can be explored in remarkable new ways. The goal of this program is to 
determine all coefficients in the superpotential, the explicit form of the Kahler potential 
and, ultimately, to perform first-principle calculations of physical quantities such as the 
relative quark and lepton masses. 

In this paper, we make substantial progress towards this goal by extending previ- 
ous work |29| to include vector bundles defined over manifolds with higher-dimensional 
Kahler cones; that is, for which /i^'^ > 1. Importantly, our results allow one to study 
arbitrary vector bundles arising in heterotic string compactifications and to determine 
whether such geometries admit N = 1 supersymmetric vacua. The problem of finding 
the Kahler cone substructure, that is, the regions in Kahler moduli space where a given 
holomorphic vector bundle is or is not slope-stable, is a notoriously difficult one. In 
particular, the difficulty of a direct stability analysis generally increases rapidly with 
the dimension h^'^ of the Kahler cone. One of the great advantages of the algorithm 
presented in this paper is that, unlike a standard analytic analysis, our numerical cal- 
culations can be performed with essentially equal ease in arbitrary /i^'^. This provides 
an important new tool in the study of supersymmetric heterotic vacua. 



The structure of the paper is as follows. To begin, in Section 2 we provide a brief 



review of the numerical algorithm for computing the Ricci-flat metric, g, and the Her- 



mitian Yang-Mills connection, A. Starting in Section 3, an overview of Donaldson's 
algorithm for computing the Ricci flat metric [20l [21j on a Calabi-Yau manifold is 
given. In particular, the numerical implementations developed in |301l31j and [32tl24j are 



discussed. Next, in Subsection 3.2, we outline the recent generalizations of Donaldson's 
algorithm presented in |231 [2^ [29] . These make it possible to compute Hermite-Einstein 
metrics on holomorphic vector bundles over a Calabi-Yau manifold and, hence, to solve 
for the unique gauge connection satisfying the conditions for A/" = 1 supersymmetry. 

Both the original Donaldson algorithm and its generalizations to connections rest 
on finding a particularly "nice" projective embedding. In the case of the Ricci-flat 
metric on a Calabi-Yau manifold, the embedding is defined from X into some higher- 
dimensional projective space via the global sections of some ample line bundle, C®^^ , on 
X. In the case of the connection on a rank n bundle V, a map into the Grassmannian 
G{n,Nku - 1) is constructed out of the Nk^ sections oi V ® C®^h ^ where C^^" is 
some ample line bundle on X. Using either one of these embedding^ a metric can be 
pulled-back to the Calabi-Yau manifold and vector bundle, respectively. One obtains a 
feg-dependent sequence of Kahler metrics X and a /sj^-dependent sequence of Hermitian 
fiber metrics on V. The degrees of freedom of every embedding parametrize a family of 
pulled-back metrics. By tuning the embedding to the so-called "balanced embedding" 
for each degree kg and kn, the Kahler metrics converge to the Ricci-flat metric on X 
and the Hermitian fiber metrics converge to a Hermite-Einstein metric on V. Finding 



^Technically, we only require a map that is an immersion rather than an embedding. However, we 
need not make the distinction in this paper. 
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the balanced embedding is solved by Donaldson's T-operator and its generalization due 
to Wang and others [IllIHalEOlEIlEaEallMlES]. Roughly, the T-operator acts on 
embeddings of fixed degree and has the balanced embedding as a fixed point. For 
a Calabi-Yau manifold, iterating the T-operator will always converge to a balanced 
embedding, and, therefore, to the Ricci-flat metric on X in the limit that kg — >• oo. In 
the case of the connection, the iteration of the T-operator for fixed k}{ is not guaranteed 
to converge. In fact, it converges to a balanced embedding if and only if the bundle V 
is Gieseker-stable. Furthermore, when V is slope-stable then the sequence of balanced 
embeddings define a fiber metric converging to the Hermite-Einstein fiber metric in the 
limit that kn — >• oo. 

In [Section "4 we modify a number of the numerical tools developed in previous 
work [29j to enable us to compare the convergence of the generalized Donaldson algo- 
rithm for different rays (or "polarizations") in Kahler moduli space. These results are 
illustrated in Subsection 4.2 by an indecomposable, rank 2 vector bundle defined over 
the K?> surface via the monad construction [36 ^ I37 t l38]. Furthermore, in Section 5, we 



present the technical details of a newly developed, rapid numerical scheme for integrating 
over a Calabi-Yau manifold. 

Section 6] provides a criterion to decide whether or not a vector bundle is slope-stable 



for a given polarization, without the need to explicitly compute the connection. Hence, 
this check can be rapidly applied to decide slope stability. In particular, we use a result 
of Wang [23j which states that the iteration of the T-operator will reach a fixed point 
if and only if the defining vector bundle is Gieseker-stable. While Gieseker stability 
is not sufficient to guarantee a solution to the Hermitian Yang-Mills equations (and, 
hence, a supersymmetric heterotic vacuum), it still provides valuable information. In 
particular, while a slope-stable bundle is automatically Gieseker stable, the converse 
does not follow. A Gieseker-stable bundle need only be slope semi-stable. Despite these 
subtleties, we extract results from the T-operator convergence which can be used to 
determine the Kahler cone substructure. Related to the question of semi-stability, we 
consider the dependence of slope-stability on the vector bundle moduli H'^iV (8) V^) in 



Section 7 In particular, the numerical algorithm is tested on a "stability wall" j39t l4Uj in 



Kahler moduli space, the boundary between slope-stable and unstable regions. We find 
that the generalized Donaldson algorithm is sensitive to the bundle moduli dependence 
and, hence, our results also distinguish the marginal cases of slope poly-stable bundles 
from strictly semi-stable ones. 

In Section 8, we extend our study to higher-dimensional spaces by presenting an 
example of Kahler cone substructure of a rank 3 monad bundle defined over a Calabi- 
Yau threefold constructed as a complete intersection in a product of projective spaces. 
We conclude and discuss future work in [Section 9[ 



2 Hermitian Yang-Mills Connections and Fiber Metrics 

A supersymmetric x heterotic string compactification is specified by 1) a com- 
plex d-dimensional Calabi-Yau manifold, X, and 2) a holomorphic vector bundle, V, 
with structure group K C Eg defined over X. The gauge connection. A, on V with 
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associated field strength, must satisfy tfie well-known Hermitian Yang-Mills (HYM) 
equations [9]. For general U{n) structure groups, these equations are given by 

Fij = F\j = 0, g'^F^ = fi{V) ■ (2.1) 

where g^^ is the Calabi-Yau metric, n is the rank of V, the scalar /xCV) is a real number 
associated with V and i, j = 1, . . . d, run over the holomorphic indices of the Calabi-Yau 
d-fold. Our primary interest is in Calabi-Yau threefolds, since compactification on these 
give rise to A/" = 1 supersymmetric theories in four dimensions. However, in order to 
present simple illustrations of the techniques introduced in this paper, we will discuss 
Calabi-Yau twofold {K3) as well as threefold examples. It is not strictly necessary for 
the first Chern class of the bundle to vanish [IT] , and the methods used in this paper 
would work just as well in that setting. However, most realistic compactifications are 
based on structure groups K = SU{n) C U{n), and these will be our main focus. When 
K = SU{n), the parameter /xCV) = and eq. (|2.1[) reduces to 



Fi, = Fij = 0, g'^F^ = 0. (2.2) 



While eq. (2.2) are the relevant equations for realistic heterotic compactifications, math- 
ematically it will often be useful to discuss the Hermitian Yang-Mills equations in full 
generality. 



A solution to (2.1) is equivalent to the bundle V carrying a particular Hermitian 
structure. An Hermitian structure (or Hermitian fiber metric), G, on V is an Hermitian 
scalar product Gx on each fiber V(x) which depends difFerentiably on x. The pair 
(V, G) is often referred to as an Hermitian vector bundle. For a given frame, ea{x), the 
Hermitian structure specifies an inner product as 

iea,eb) = Gab, G = Gl (2.3) 

A choice of frame provides the necessary coordinates to express the covariant derivative 
in terms of the connection, 

Div'^ea) = {dv^)ea + v^A\eb. (2.4) 

Imposing compatibility of the connection with the holomorphic structure of the bundle 
and the fiber metric determines the connection uniquely up to gauge transformations. 
Written in the most useful gauge choice for our purposes, the connection is 

^ = 0, A = G-^dG. (2.5) 

One can then rephrase the Hermitian Yang-Mills equation for F^^'^^ in (2.1) as a condi- 



tion on the bundle metric, 

/u(V) • l„xn = /'F,- = g'^d-jA, = g^%{G-^diG). (2.6) 

A metric G on the fiber of V satisfying this equation is called an "Hermite-Einstein 
metric" . By integration, this metric can be used to define an inner product on the space 
of global sections of V, s% where a = 1, . . . hP[X^ V), 



sp) = I slGb-asl dVol. (2.7) 

X 
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The above notions in differential geometry can be related to seemingly very different 
concepts in the algebraic geometry of holomorphic vector bundles. Relating the two 
approaches has made it possible to better understand both. For Kahler manifolds, the 
relationship can be summarized as follows: 

Theorem 1 (Donaldson-Uhlenbeck-Yau |1H I12j). On each slope poly-stable holomor- 
phic vector bundle, V, there exists a unique connection satisfying the general Hermitian 



Yang-Mills equations eq. (2.1). Moreover, such a connection exists if and only if V is 
slope poly-stable. 

Thus, in the heterotic string context, to verify that a gauge vector bundle is consistent 
with supersymmetry one need only verify that it is slope poly-stable. The notion of 
slope-stability of a bundle J-' over a Kahler manifold X is defined by means of a real 



number (the same which appeared in eq. (2.1)), called the slope: 



where d is the complex dimension of the Kahler manifold. Here, uj is the Kahler form 
on X, while rk(J^) and ci(J^) are the rank and the first Chern class of J- respectively. 
A bundle V is called stable {semi-stable) if, for all sub-sheaves T CV with < rk(J^) < 
rk("\7), the slope satisfies 

^(V). (2.9) 

A bundle is poly-stable if it can be decomposed into a direct sum of stable bundles which 
all have the same slope. That is, 

v = 0v„, ^l{v,) = ^i{v). (2.10) 

n 

From the above definitions, it is clear that the condition of slope-stability on a Calabi- 
Yau manifold depends on all moduli of the heterotic compactification. To be specific, 
consider a Calabi-Yau threefold. Here, the moduli are the h^'^{X) Kahler moduli, the 
/i^'^(X) complex structure moduli, and the /i^(End('V)) vector bundle moduli. The 



dependence on Kahler moduli is explicit in eqns. (2.8) and (2.1). Since slope stability 
is an open property 02], it depends only on a Kahler form, co, defined up to an overall 
scale. We refer to this one-parameter family of Kahler forms (which define a ray in 
Kahler moduli space) as a choice of "polarization" and frequently make no distinction 
between a particular uj and its associated polarization. It is possible to expand the 



Kahler form uj in (2.8) as a; = f^ujr, where uir are a basis of (1, l)-forms and f are the 
real parts of the Kahler moduli. Written in terms of the triple intersection numbers drst 
of the threefold, the slope is simply 

^(V) = ^ E drsMVytHK (2.11) 

^ ^ r,s,t=l 

The complex structure moduli of the Calabi-Yau manifold and the vector bundle moduli 
enter through the notion of a subsheaf C V. Thus, finding a solution to the Hermitian 
Yang-Mills equations, or determining whether the bundle is slope-stable, is a question 
that must be asked after selecting a particular point in moduli space. 
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3 The Generalized Donaldson Algorithm 



Many of the challenges associated with string compactifications on a Calabi-Yau d- 
fold X arise from the difficulty in determining the explicit geometry. The simplest 
J\f = 1 supersymmetric vacuum solutions require a Ricci-flat metric, i^jj, on X and a 



Hermite-Einstein bundle metric, Gab, satisfying (2.6) as discussed above. While Yau's 
theorem [lUj ensures that a Ricci flat metric exists on a Calabi-Yau manifold, and the 
Donaldson-Uhlenbeck-Yau theorem [12] provides for the existence of a Hermite- 
Einstein metric on a slope-stable bundle, no analytic solutions for either the metric or 
connection have yet been found. 

However, recent work has made it possible to find accurate numerical solutions for 
both metrics and connections. An algorithm was initially proposed by Donaldson for 
the computation of Ricci-flat metrics [191 ISQl |21] , and was implemented numerically and 
extended in[l3l[3nil3ll|32l[2ll[271[29]. What we refer to as the "generalized Donald- 
son algorithm" is an extension of Donaldson's approximation scheme which numerically 



approaches an Hermite-Einstein bundle metric, solving (2.6). This was developed math- 
ematically in |22[ 123] and implemented numerically in [24, 29j. A thorough review of 
the Donaldson algorithm and its extensions is beyond the scope of this paper. We re- 
fer the reader to [29j for more details. However, in order to proceed with our present 
investigation of Kahler cone substructure, we provide here a brief review of the central 
ingredients of the (generalized) Donaldson's algorithm and set the notation that will be 
used throughout this work. 

3.1 Donaldson's Algorithm 

We begin with an overview of Donaldson's algorithm for approximating the Ricci flat 
metric on a Calabi-Yau manifold. The first ingredient we need is one of the simplest 
Kahler metrics, the Fubini-Study metric on F"". This is given by Qpsij — ^didjKps, 
where 

KFs = -lnS2h'~^ZiZj (3.1) 

and is any Hermitian, positive, non-singular matrix. 

Since it is always possible to embed X C P" for some large enough n, the Fubini- 
Study metric can be used to induce some metric on any Calabi-Yau manifold X. Such 
a metric will not be Ricci-flat, for otherwise one could easily write down an analytic 
expression for the Calabi-Yau metric. It is tempting to wonder whether there exists a 



generalized version of eq. (3.1) with enough free parameters to provide a more versatile 
induced metric on X? The central idea of Donaldson's algorithm is to find such a 
generalization and a procedure for successively tuning its free parameters to approximate 



the Ricci-flat metric. The obvious generalization of eq. (3.1 ) is to replace the degree one 
polynomials with polynomials of higher degree. That is, 

K = i- In _ f^'-'^^'-'^z.. . . ■ ... % (3.2) 

ji...ifeji...jfc 
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where h^i---^kji---jk Hermitian. This new Kahler potential now has (n + l)^'^ real pa- 
rameters. This generalization can, in fact, be seen in a more systematic way by using 
holomorphic line bundles over X. The Kodaira Embedding Theorem |44j tells us that 
given an ample holomorphic line bundle L over X with ni = h^{X,L) global sections, 
one can define an embedding of X into projective space via the sections of C'' = 
for some k. That is, choosing a basis for the space of sections, Sa G H^{X,L^) where 
< a < nfc — 1, allows one to define a map from X to W'''~^ given by 

ifc: X^P"^-i, (xo,...,Xrf_i)^ [so(x) : ... (3.3) 

where Xi are holomorphic coordinates on the Calabi-Yau manifold. If C is sufficiently 
ample, eq. (3.3) will define an embedding of X C P"*^"^ for all with k > ko for some 
ko- 

In terms of this embedding via a line bundle C, one can view the generalized Kahler 
potential in eq. (3.2), restricted to X, as simply 

ni^ — l 



fevr 

a,l3=0 



Geometrically, (3.4) defines an Hermitian fiber metric on the line bundle £ itself. It 
provides a natural inner product on the space of global sections 



where 

dVolcY = A n (3.6) 

and $7 is the holomorphic (3,0) volume form on X . 

With the initial Kahler metric eq. (3.4) in hand, we must now proceed to system- 
atically adjust it towards Ricci fiatness. To accomplish this, the notion of a balanced 
metric is required. Note that, in general, the matrices h°'^ and in (3.5) are com- 
pletely unrelated. However, for special metrics, they may coincide. The metric h on the 
line bundle L is called balanced if 

(M,^)-i = h"^. (3.7) 

Donaldson first recognized that balanced metrics lead to special curvature properties. 
These can be summarized as follows [191 1201 [2H [M]: 

Theorem 2 (Donaldson, Keller). For each k > 1, the balanced metric h on exists 
and is unique. As k ^ oo, the sequence of metrics 

a,fS=0 

on X converges to the unique Ricci-flat metric for the given Kahler class and complex 
structure. 
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The central task of Donaldson's algorithm is thus to find the balanced metric for 
each k. To this end, Donaldson defined the T-operator as 



For a given metric h, it computes a matrix T{h). If this matrix equals M^j^, we have a 
balanced embedding. To find this fixed point, simply iterate (3.9) as follows. 



Theorem 3. For any initial metric ho (and basis Sa of global sections of L^), the 
sequence 

hm+l = {T{hm)y^ (3.10) 

converges to the balanced metric as m —)■ oo. 

Happily, in practice, very few (~ 10) iterations are needed to approach the fixed point. 



(k) 

Henceforth, we will also refer to g:^ in eq. (3.8), the approximating metric for fixed k, 



as the balanced metric. It should be noted that, to find the balanced metric at each 

we 



step k, one must be able to integrate over the Calabi-Yau threefold. In Section 5 
will discuss the new adaptive mesh numerical integration scheme used throughout this 
work. 

As one final ingredient in the algorithm, one must be able to quantify how closely the 
numerical metric approximates the Ricci-flat metric. A variety of such error measures 
were given in [29]. Recall that, given an sufficiently ample line bundle C, one can find 
a Kahler form 

= ^g^^^dzi Adzj (3.11) 

corresponding to the balanced metric associated with C^. Note that the Kahler class of 
this Kahler form is [u!k] = 2ttci{C^) and the associated volume is 

1 /■ 4 



Vol, = jj (3.12) 



X 



where denotes the (d, d) volume form w A • • • A w. 

In this paper, we will measure convergence of the Donaldson algorithm via the Ricci 
scalar in 

ll^i^lU = Vof-'^^/^y \Rk\ Vdet^fe d^'^x. (3.13) 

On a Calabi-Yau manifold, H-B-ffUfc = 0{k^^) as A; — )• c« and, hence, this error measure 
should approach zero. As a final note, we will henceforth denote the degree of twisting, 
given by the integer k in C^, as kg to make it clear that this integer is associated with 
the computation of the metric. 



A summary of Donaldson's algorithm for Ricci fiat metrics is provided in Table 1 
We now turn to the generalized Donaldson algorithm for computing Hermite-Einstein 
fiber metrics on holomorphic vector bundles. 
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3.2 Hermite-Einstein Bundle Metrics 



As we saw in the previous section, Donaldson's algorithm is a powerful tool for numer- 
ically approximating the Ricci-flat metric on a Calabi-Yau manifold. In this section, 
we investigate a generalization of these techniques which can be used to approximate 
the field strength F^^'^^ of a holomorphic connection which satisfies (2.1). As discussed 



m 



Subsection 3.1[ Donaldson's algorithm for Calabi-Yau metrics can be viewed as a 



method for numerically obtaining a particular Hermitian structure on the ample line 
bundle This balanced fiber metric on allows one to define a balanced embedding 
of the Calabi-Yau space X into P"fc^i. By mapping the coordinates x £ X into the 
global sections Sa S H^{X,C^), that is. 



X h 



[so{x) : ■■■ : 



(3.14) 



we produced a map ik '■ 

X p«fe-i where Uk = /i°(X,£^). The pull-back of the 
associated Fubini-Study metric was shown in Subsection 3.1| to converge to the Ricci-fiat 
metric on X in the limit that k ^ oo. Viewed in terms of Hermitian fiber metrics on line 
bundles, it is a natural question to ask whether Donaldson's algorithm could be extended 
to develop an analogous approximation to Hermitian metrics on higher rank vector 
bundles. In particular, could one find an approximation scheme to produce an Hermitian 
metric on an arbitrary stable bundle V of rank n such that it satisfies condition (2.6)? 



Fortunately, precisely this question has been addressed in the mathematics literature |23J 
and implemented for physics in |24l [29] . 

To generalize Donaldson's algorithm, consider defining an embedding via the global 
sections of a twist of some holomorphic vector bundle V with non-Abelian structure 
group. That is, consider a map 



X h 



(3.15) 



from X £ X into the global sections S'^ G H^{X,'V ® C^), where a = ... A^^ — 1 indexes 
the hP{X,'V®C'^) global sections and the index a = 1, ... n is valued in the fundamental 
representation of structure group K (1 U (n) of the rank n bundle V. We hope then to 
define the embedding 

X ^G{n,Nk-l), (3.16) 

where G{n, — 1) denotes the Grassmannian of the relevant dimensiorj^ By the 
Kodaira embedding theorem j35|, given a holomorphic vector bundle, V, and an ample 
line bundle, C, there must exist a finite integer ko such that, for any k > ko, the twisted 
bundle V{k) = V ^ defines an embedding, ik : X ^ G{n, Nk — I). 

As in the Abelian case in the previous section, one can attempt to use this embedding 
to define a Hermite-Einstein bundle metric on V ^ C'^ and, hence, an Hermitian Yang- 



Mills connection as in (2.5) and (2.6). If £ is ample then, for some sufficiently large k, 



In this language, the Abehan case in (3.31 is simply an embedding X — >■ G(l,nfc — 1). 
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V ® will be generated by its global sections. That is, it will define an embedding as 



in (3.16). In our search for a solution to the Hermitian Yang-Mills equation (2.1), the 
connection on the twisted bundle 'V(^C^ will be closely related to the original connection, 
since such a twist only modifies the trace part of the field strength. Stated in terms of 
algebraic geometry, the process of twisting will not modify the slope-stability properties 
of "V since V ® is stable if and only if V is. 



As at the beginning of [Subsection 3.H where we chose the trial form of the Kahler 
potential in 



G in eq. (2.3 



3.4), here we begin with another simple anzatz for the Hermitian structure 
). Consider the matrix 



-l\ab 



E 

Q,/9=0 



(3.17) 



where H°^^ is a Hermitian matrix of constants and are the global sections ofV ^ C^. 
As in ( |2.7| ), this fiber metric induces an inner product on the space of sections H^{X, V0 
C^) = spanjS'a} via 



/ Q I Q \ 



Volcy 



S^{G''')-'sU\olcY 



X 



Volcy 



X 



(3.18) 

With this definition of the inner product on sections, one can give a natural generaliza- 



tion of the T-operator eq. (3.9). This generalization, 



Volcy 



Sa(s^Hs) Sn dVolcy, 



(3.19) 



X 



was introduced in [23j and studied numerically in [2 



5] . Note that if "V is a line bundle 
then eq. (3.19) reduces to (3.9) and one recovers the case of a balanced embedding 
into p^fc^i. As in the previous section, we will now describe how the iteration of the 
generalized T-operator can produce a fixed point which describes an Hermite-Einstein 
bundle metric. 

To do this however, we must introduce one additional notion of stability, namely 
that of "Gieseker stability" [42] . Let C be an ample line bundle and be a torsion- free 
sheaf. The Hilbert polynomial of with respect to C is defined as 



rk(^) 



(3.20) 



where ® C^) is the index of (8) Given two polynomials / and g, we will write 
f -< g a f{n) < g{n) for all n ^ 0. Then a bundle V is said to be Gieseker stable if, for 
every non-zero torsion free subsheaf cV, 



With this definition in hand, it was shown in [23j that 



(3.21) 
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Theorem 4 (Wang). A bundle V is Gieseker stable if and only if the k-th embedding, 
defined by V as in eq. (3.16), can be moved to a "balanced" place. That is, if there 



exists an orthonormal section-wise metric on the twisted bundle such that 

'1 TTai3 



{T{H) 



is a fixed point of the generalized T- operator . 



(3.22) 



We can use this special metric on V ^ jC'^ to define an Hermitian metric on V itself. 
Let Gc denote the balanced metric on £., and G^'^^ the balanced metric on V(8)£'^. Then 

Gk = G^^^ ^ Gl^ (3.23) 

is an Hermitian metric on V. This appears in the following important theorem 



Theorem 5 (Seyyedali, Wang). Suppose V is a Gieseker stable bundle of rank n. If 
Gk^G ^ as k — y oo, then the metric G^o solves the "weak Hermite-Einstein equation" 



R-R 



(3.24) 



where 



• R is the scalar curvature. 

• R = f R^/deUj d^'^x is the averaged scalar curvature, which is zero for any Kdhler 
metric on a manifold of vanishing first Chern class. 

We will, henceforth, denote the degree k of the embedding defined above as kjj, to 



make clear its association with the Hermitian matrix in eq. (3.17) and distinguish it 
from kg. Procedurally, the process of obtaining the Hermite-Einstein fiber metric on 
a slope-stable bundle V is very similar to that outlined for the Ricci-fiat connection 
in Subsection 3.1: for each value kn of the twisting, we iterate the T-operator asso- 
ciated with the embedding defined by H^{X,'V (g) C^") until a fixed point is reached. 



Then, by Theorem 5, the induced connection approximates solutions to eq. (3.24) as 
kn — >• oo. However, there is an immediate and important difference between this gener- 
alized algorithm and Donaldson's algorithm for Ricci-fiat metrics. While all Calabi-Yau 
manifolds admit a Ricci-fiat metric, not all holomorphic vector bundles will admit an 



Hermite-Einstein metric satisfying (2.6). That is, if one applies the algorithm to a bun- 



dle that is not slope-stable, it will not converge to a solution of the Hermitian Yang-Mills 



equations, (2.1). Moreover, it should be noted that while all slope-stable bundles are 
the converse does not hold: 



Gieseker stable ^], the converse does not hold: not all Gieseker-stable bundles are 
slope-stable. That is, there exist cases where the iteration of the T-operator does con- 
verge for fixed kn, but the sequence of metrics does not converge towards a solution of 
the Hermite-Einstein fiber metric. However, if V is a slope-stable holomorphic bundle, 
then the iteration Hm+i = T[Hm)~^ will converge at each kff, and in the limit that 
kn — >• oo, produce the Hermitian bundle metric Goo satisfying (3.24) via its associated 
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field strength defined in (2.5) and (2.6). Moreover, in the case where the Calabi-Yau 
metric g^^ is Ricci-flat, (3.24) simply reduces to (2.1). Thus, we have found a solution 



to the Hermitian Yang-Mills equations. 

However, one must be careful. Despite having found a Hermite-Einstein bundle 
metric (and, hence, HYM connection) associated with the twisted bundle V C]^" , and 
an Hermitian metric Goo satisfying (3.24), our task is not yet complete. We still need 



to explicitly determine the connection on the bundle V itself satisfying the Hermitian 
Yang-Mills equations, (2.1). Since the process of twisting V by a line bundle CJ^" in the 



above construction clearly modifies the trace-part of the connection, one must subtract 
this line bundle contribution to get the connection on V only. To do this, we have to 
separately find a suitable metric Gc on £. For example, one could compute the balanced 
metric '"^ = s^hs on for some sufficiently large k^. Then Gn = [s^hsY^^'^ would 



approximate the constant curvature Hermitian fiber metric on L and, as in eq. (3.23), 
we find that 



G 



( S^HS] ( s^hs 



(3.25) 



is the fiber metric ([2^ on V. As before, S G H^{X,V (g) L^^) and s G H°{X,C'"') 



are the relevant global sections. Using eqns. (2.5) and (2.6), in terms of the Hermitian 
metric, the connection on V is then given by 



A{V) = d 



S^HS] ( s^h 



\ k„/kh 



kh 



(3.26) 



That is, one can "untwist" the connection simply by subtracting the trace of the Abelian 
connection on C''^ to produce the U{n) connection on V. The curvature is given by 



^(0,2) ^ ^(2,0) ^ 0^ ^ g'~^djAi = g'~^djdi In (s^Hs) (s^h 



-kn/kh 



(3.27) 



As shown in |29], when V is a U{n) bundle, the most efficient way to perform this 
untwisting is not by computing an independent balanced metric (Gc)^'^, but directly 
using the induced Hermitian fiber metric on the determinant line bundle A'^{V (>S) L^^) 
of V C^" . In particular, we choose kh = Ta,'nk{y)k}j . It follows that the Hermitian 
metric on C^*^ is 



{Gcf^ = det(G('=«)) = det {S^HS). 



(3.28) 



Let A'-'^^^ be the eigenvalues of g^^ F^-^"^ on V , and let A be the corresponding 



eigenvalues of g^-' Ff: on V after untwisting. Using eq. (3.27), we obtain that 



A,; = A 



rank V *J 



{kH) 



A, 



{kH) 



rank V 



(3.29) 



where i = 1, . . .n where n is the dimension of the fundamental representation of the 
structure group of "V. Therefore, the effect of this untwisting is precisely to subtract, at 
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each point, the average of the eigenvalues. Hence, in 
as subtracting the trace. 



29] we referred to this untwisting 



The eigenvalues in (3.29) are a pointwise measure of the error in the numerically 
derived connection. For a slope-stable bundle, Aj — >• ^ as one increases kn ^ oo. To 
properly define an error measure for the approximation to the Hermitian Yang-Mills 



connection, we must test the approximation at all points and, hence, integrate (3.29) 
over X. As in [29], we define the error measure 



27r Volfc„ rank(V) 



X 



(E 



(3.30) 



where Vob 



is the volume computed in (3.12). For a slope-stable bundle, T{Ay) — )• 
oo. This is simply a global check of the eigenvalues in (3.29). To 



rank(V)/u as kn 

summarize the results of this section, the generalized Donaldson algorithm for numeri- 



cally approximating a Hermitian Yang-Mills connection is presented in Table 1 



4 Kaliler Cone Substructure 

4.1 Modifications For Higher Dimensional Kahler Cones 

One of our central motivations in this work is to understand the generalized Donaldson 
algorithm on manifolds with higher dimensional Kahler cones, /C, that is, dim(Pic(X)) > 
1. In particular, we will compare the behavior of bundles under the algorithm for differ- 
ent choices of polarization. In general, holomorphic vector bundles can display different 
slope-stability properties for different choices of polarization, that is, along different 
rays in the Kahler cone jf] That is, a given bundle may be slope-stable in some sub- 
cone /Cstabie C /C, but be slope-unstable (and, hence, break supersymmetry) in other 
sub- regions /Cunstabie C /C. This substructure is of interest both mathematically and 
physically, with applications ranging from supersymmetry breaking in heterotic J\f = 1 
supersymmetric vacua [ 4:6\ \39\ BOl HTJ HH] H9] to the computation of Donaldson-Thomas 
invariants on Calabi-Yau threefolds \50 \ [51 } [52]. In general, it is a difficult task to deter- 
mine the global slope-stability properties of a vector bundle, V, throughout the Kahler 
cone. In particular, this analysis scales badly with the dimension of Kahler moduli space, 
h^'^{X). Already for h^'^ > 4 it becomes prohibitively difficult to analytically analyze 
the stability of a bundle except in special cases. As a result, it is of considerable interest 
to ask the question: Can the generalized Donaldson algorithm provide an efficient probe 
of Kahler cone substructure and vector bundle stability for higher dimensional Kahler 



cones? In principle, the connection algorithm reviewed in Subsection 3.2 shows no dif- 



ference in computational difficulty for any dimension h^'^. That is, it depends only 



on a one-dimensional ray (defined by the line bundle C in Step 1 of Table 1 and the 



embedding ( 3.15| )) and not on the dimension of the Kahler cone containing that ray. As 



we will see in the following sections, the generalized Donaldson algorithm does indeed 
provide a powerful new tool for analyzing Kahler cone substructure. 



^The radial direction along a fixed ray only parametrizes the overall volume and does not change the 
stability properties. 
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In order to pursue this goal, however, one will need a way to compare the convergence 
of the algorithms (for both metric and connection) for different rays in Kahler moduli 
space. A number of properties change in the case that h^'^ > 1 and, in particular, a few of 
the definitions introduced in Section 2, and in the previous literature [241 '29j, need some 
modifications in order to make sensible comparisons for different polarizations. One of 
the first of these is the way in which we measure the complexity of the embeddings in 
(3.3) and (3.15). Recall that the algorithms described in Section 2 rely on defining an 
embedding into some high-dimensional Grassmanian, (3.15). For example, to compute 
the Ricci-flat metric of X, we define the embedding X — )■ P" via the global sections 
H^{X, C^s). For a manifold with h}'^ = 1, it is clear that as we increase the degree, 
kg, of twisting, we increase the number of global sections and, hence, as described in 



Subsection 3.1, the accuracy of the metric approximation. For example, in |29j we 
computed Ricci-flat metrics on the Quintic hypersurface in P^, where the global sections 
of the embedding line bundle, C = 0(1), increase with kg as H^{X, 0{kg)) = l{bkg+kl). 

However, for manifolds with /i^'^ > 1 the situation becomes more subtle if one wants 
to compare results for two different polarizations, defined by line bundles Ci and £2- 
As an example, consider the Calabi-Yau 3-fold X defined as a (2, 4) hypersurface in 
P^ X P'^. This manifold has /i^'^ = 2 and its Picard group is spanned by the restriction 
of the respective hyperplanes of P^ and P^ to X (respectively the line bundles 0(1,0) 
and 0(0,1)). Now, consider two distinct polarizations defined by Ci = 0(2,1) and 
C2 = 0(1, 2). We can define an embedding of X into some projective space using either 
of these ample line bundles. However, the sections of each grow very differently in 
H^{X,Ci^^) where i = 1, 2. These sections grow with kg as 

H\X,0{2,lf''^) = hg{23 + 13kg^), H°{X,0{l,2f'^) = hg{28 + 32kg^). (4.1) 
o o 

Hence, if we computed the metric for each of these polarizations to the same degree, 
say kg = 10, we would have very different results. From Ci we would have defined an 
embedding with 4, 410 sections, while with £2 we would have 107, 600 sections - and, 
via the algorithm of Subsection 3.1 a far more accurate approximation to the Ricci-flat 
metric. 

In order to sensibly compare results for different polarizations then, instead of the 
degree of twisting kg (or kfj in the case of the connection), we specify the number 
of sections in H^{X, C^^). When we compare the results for different polarizations, 
we compare them at orders chosen so that there is an approximately equal number of 
sections. This procedure was first introduced for Ricci fiat metrics in \30\ I31j. 

There is one further modification one must make to the definitions of ISection 21 This 
arises in error measures used to test the accuracy of the approximation to an Hermitian 
Yang-Mills connection. In order to determine whether or not a bundle is slope-stable 
for different polarizations, we must charge the normalization of our error measure in 
eq. (3.30). For example, take the case of a general rank n stable vector bundle with 
structure group U{n) satisfying the general Hermitian Yang-Mills equations (2.1) on a 
d-dimensional Kahler manifold X. Since g^^'F^j = f^^V) ■ Inxn, it is straightforward to 
see that in this case 

IX 



r(^v)= / ci(V®£'=«) Aw"'-^ gZ. (4.2) 
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Figure 1: A plot of the normalized error measure, (4.3) 
for the sum of line bundles, 0(0, 1) 0(0, —1) 
(2, 3) hypersurface in X P3. Using the error measure introduced in 



The results shown are 
on the K3 defined as a 



(3.30), the results would vary for different rays in Kdhler moduli space. 



However, with the new normalization in (4.3) the results are uniform. 



As a result, r manifestly depends not only on the first Chern class of V, but also on 
the choice of polarization oj. In fact, it jumps by an integral amount if one changes the 
polarization. In order to compare the r error measure for different choices of polarization, 
we will introduce a new normalization that will remove the polarization dependence and 
make the initial values of r more uniform for different twists V C]^^ . 

As an example, take a sum of line bundles Ca of different first Chern class 
and let fisum = ©a/^(^a)- given polarization C = 0{tr) = 0{ti,t2, ■ ■ ■) with 

r = 1, . . . , h^'^ , define Upoi = /^(/^)- We then introduce the new normalized error measure 

T{tr) = -^T{Avh,. (4.3) 
IJ'Sum 

This choice is made so that r(fj.) will be independent of the polarization and equal to 1 



for any sum of line bundles. This normalization is shown in Figure 1 for the sum of line 
bundles 

0(0,1)0 0(0,-1) (4.4) 

on the K3 surface defined via a degree (2,3) hypersurface in x F'^. Without this 
normalization the error measure of the sum of line bundles would vary according to the 
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polarization, making it hard to compare the different directions in the Kahler moduli 



space. However, as is illustrated by the figure, the error measure in (4.3) produces 
uniform results. With these new definitions in hand, we turn to our first systematic 
study of Kahler cone substructure. 

4.2 A Bundle On An Elliptic K3 Surface 

In this section, we explore the Kahler cone substructure described above in an explicit 
example. In particular, we consider the elliptic K3 surface X defined as a degree (2, 3) 
hypersurface in x P^. This representation of the K3 has dim(Pic(X)) = 2, that 
is, a 2-dimensional algebraic Kahler cone, /C. Expanding the Kahler form on X in 
a harmonic basis of (1, l)-forms as w = tiui + t2Uj2, the Kahler cone is the positive 
quadrant, ti,t2 > in the coefficients tr- Over this space, we will denote line bundles 
by 0{ti,t2), where 0(1, 0) and 0(0, 1) are the pull-back of the hyperplane bundles of P^ 
and P^, respectively. We will consider below the stability properties of a sample SU{2) 
bundle over different polarizations in the Kahler cone, /C. 

We compute both the Ricci-flat metric on the base manifold and the Hermitian Yang- 
Mills connection on a bundle for a number of different polarizations. In the case of the 
metric, this amounts to computing the fiber metric on the line bundle C = 0{ti,t2) for 
different positive choices of ti,t2- As discussed above, in order to compare the accuracy 
of these approximations to a Ricci-flat metric for different polarizations, t2/ti, one can 
no longer simply specify a twisting £^^9. Instead, as we vary the coefficients tr, we will 
compare the number of global sections that are generated by C^^ for different choices of 
C. That is, for each choice of C, we will compute the metric up to the largest degree kg 
such that there are < 500 sections in H^{X, C^^). In the calculation. 



m 



Subsection 3.1 



the metric algorithm utilized 202, 800 points in the adaptive numeric integration (as will 



be described further in Section 5) and the metric T-operator was iterated 30 times. 

On this Calabi-Yau twofold we now define a rank 2, holomorphic vector bundle with 
structure group SU{2). This sample bundle is defined through the so-called monad 
construction pj [371 l3Hl l53] . 

0^O(-2,-l) ^ O(-2,0)®2©O(2,-l) ^V^O. (4.5) 

Here V is defined as the cokernel of a generic map / with bi-degrees ((0, 1), (0, 1), (4, 0)) 
between the direct sums of line bundles. Using the techniques of |421 154j. it is straight- 



forward to prove that V in (4.5) is destabilized in part of the Kahler cone by the rank 
1 sheaf 

^-F^ O(2,0)®2 ^ 0(2,1) ^ (4.6) 

with ci{J-) = (2, —1). The intersection numbers drs of the two hyperplane classes in the 
(2, 3) K3 surface are du = d2i = 3, d22 = 2, and du = 0. Using the definition of slope 



in (2.8), it can be verified that V is slope-stable when ti/t2 > 4/3 and unstable when 



ti/t2 < 4/3. That is, the Kahler cone exhibits the substructure 

(4.7) 



^stable = {0{tl,t2)\ ^ > |}, 



/Cunstable = {0(*l,i2)| fr < §}• 
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Figure 2: The Kdhler cone substructure associated with the SU(2) bundle in (4.51 



The normalized error measure, (4.3), is shown for different choices of 
polarization in the Kdhler cone. The presence of V clearly divides the 
Kdhler cone, JC, into two chambers, ICstabie o-nd ICunstahie, corresponding 



to the stable/unstable regions described in (4.7) 



Hence, in the stable region we expect the T-operator to converge, whereas in the unstable 
region we do not expect to find a fixed point. 

As discussed in Subsection 3.2, to apply the generalized Donaldson algorithm one 
must define the embedding : X ^ G{n, Nf^^j). To do this, one must compute the 
global sections of the twisted line bundle V (8) C'^^ for some ample line bundle £. For 
the bundle defined in (4.5), the global sections H^{X,'V (g) C''") can be computed for 
any choice of twisting. Multiplying eq. (4.5) by C'^" = 0{ti,t2)^" , we obtain the short 
exact sequence 



1) 



^ OikHh - 2, kHt2 - 1) Oiknti - 2, fcHt2)^ 0{kHti + 2, kHt2 - 
Then the global sections are given simply as the cokernel 

H\X, OikHh - 2, kHt2)®^ e OikHh + 2, kHt2 - 1)) 



H^iX,V(S)Oih,t2)''"] 



fimix,oikHh-2,kHt2-i))) 



0. (4.8) 



(4.9) 



where both parts of this quotient are the global sections of sums of ample line bundles 
when knh > 2 and kHt2 > 1. 
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For V in (4.5), we will apply the generalized Donaldson algorithm as reviewed in 

In particular, we will compute Hermitian bundle metrics 



Subsection 3.2 and Table 1 



on V and determine whether they converge to the Hermitian Yang-Mills connection 
for different choices of polarization. Specifically, we will compare the results for six 
different choices of polarization: three lying within /Cstabie; two in /Cunstabie and one on 
the boundary defined by ti/t2 = 4/3. The results are shown in Figure 2, where we plot 
the normalized error measure, eq. (4.3), for the Yang-Mills connection in various 
directions in the Kahler cone. 

32 



= O(-2,0)®2eO(2,-l) 



Since V is a cokernel associated with the direct sum 
the normalization fisum in eq. (4.3) was chosen to be 
IJ-sum = 2/i(C'(— 2, 0)) + fj,{0(2, —1)) in order to meaningfully compare the different rays 
in the Kahler cone. In the stable sub-cone, the eigenvalues of g^'' F^j can clearly be seen 
to be approaching zero as one increases the number of sections in (X, V ® C^" ) , as 
expected. In the unstable region there is no such convergence. Both observations are in 
perfect agreement with eq. (4.7). The connection on V was computed with Nq = 74, 892 
points (adaptive) and the connection T-operator was computed with 100 iterations at 



each graph point in Figure 2 



The results of this section clearly indicate that the generalized Donaldson algorithm 
can be used to investigate Kahler cone sub-structure. We will explore this in more detail 
in the following sections, but first we will provide a description of the novel integration 
method implemented and used throughout this work. This integration scheme provides 
a significant increase in computation speed and makes it possible for us to analyze a 
wider range of examples. 



5 Adaptive Integration 

5.1 Rectangle Method vs. Monte-Carlo 

All approaches to numerical geometry of Calabi-Yau threefolds, be it Donaldson's algo- 
rithm tia[2Ql[2llll3l[30l[3ll[32l[23[271[29],it^ generahzation to Hermitian Yang-Mills 
bundles or direct minimization [27] all use a spectral representation of the 

geometric data. That is, the tensors describing the geometry are eventually expanded 
in a suitable basis of functions, and the problem reduces to finding the "best fit" coeffi- 
cients. This is in contrast to the traditional finite elements methods, where one directly 
discretizes spacetime. As a general rule, finite elements work well in low dimensions, 
but spectral representations are necessary in high-dimensional problems. 

The basic numerical step that every spectral algorithm relies on at its core is to 
integrate over the base manifold. This is necessarily so because only by evaluating the 
spectral basis everywhere on the manifold can one draw conclusions about the global be- 
havior of the geometric object of interest. The most straightforward integration scheme 
is to split the integration domain into equal-sized pieces and use the multidimensional 
generalization of the rectangle rule. In practice, the volume elements can only be chosen 
of equal size with respect to an auxiliary metric. But this is then easily corrected for by 
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weighting the individual points accordingly, 



^d'^x = w{x) y/^d'^x w{x) = ^0=. (5.1) 

y 5aux 

Here, the scalar weight function w{x) is simply determined by the auxiliary measure 
(given by the point distribution) and the desired measure d'^x. The Calabi-Yau case 
is particularly simple, since the Calabi-Yau volume form ^gcY d'^x = Q Aft is known 
analytically. 

The disadvantage of this direct approach is that it requires extensive knowledge about 
the geometry and topology of the manifold to construct a constant (auxiliary) volume 
cell decomposition. This is why it has only been used for complex surfaces [551 l56] . 
that is, real 4-dimensional manifolds. A solution to this problem was devised in |32| . 
where a Monte-Carlo integration scheme was proposed that needs as input only the 
defining equation of a Calabi-Yau hypersurface (or complete intersection). The key to 
this approach is that one can generate random points with known distributions using 
zeroes of random sections of line bundles [571 EH]- As an example, consider the quintic 
i : Q — )• embedded in projective space. A line C is defined by three linear equa- 
tions in the homogeneous coordinates, that is, three sections of 0p4(l). Any one linear 
equation is defined by its 5 coefficients, so one can talk about random sections with an 
S'[/(5)-invariant distribution. Three linear equations intersect the quintic hypersurface 
in 5 points, and by the general theory the probability distribution has measure i*(a;Fs)^- 
The Calabi-Yau volume form is given by the residue integral 

in a local (holomorphic) coordinate patch {pi, P2, P3, Pi), thus determining the weight 
function w{x) for the random points. The disadvantage of this Monte-Carlo integration 
scheme is that one is forced to use the random point set without modification. In partic- 
ular, it tuns out that the point weights w{x) fiuctuate over a large range with the error 
accumulating in the badly-sampled regions. This was mitigated in [59] using stratified 
sampling, at the cost of having to work with higher-dimensional Kodaira embeddings. 



5.2 An Adaptive Integration Algorithm 

For the purposes of this paper, we developed a combination of the best features of the 
naive higher-dimensional rectangle rule and the Monte-Carlo sampling. In a nutshell, 
the idea is to first parametrize point (s) of the Calabi-Yau threefold similarly to the 
parametrization by sections of three line bundles. Then integrate using the standard 
higher- dimensional rectangle rule by constructing a suitable cell decomposition of the 
parameter space, not the Calabi-Yau manifold. 

As the simplest example, consider the (2, 3)-hypersurface in X C P^ x P^, which is 
a K3 surface. The projection vr : X — )• P^ onto the P^ factor is, generically, two-to-one. 
Locally, there are two inverses vr^^, vTg"^ for the two-sheeted cover tt. With it, one can 
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rewrite the integration as 



2 



dz 



, V^aux d z, (5.3) 

V ffaux 



Wi{z) 



where the sum runs over the different sheets. The weights Wi{z) are the product of 
the Jacobian of the coordinate transformation and, as before, the scalar factor required 
to transform the auxiliary measure into the desired Calabi-Yau measure. Finally, it 
is easy to integrate over projective space. For our purposes, we will use an adaptive 
integration scheme where we start with a decomposition of into cells of equal volume 
with respect to some convenient auxiliary measure. Then, if any of the weights Wi{z) is 
significantly larger than the average weight, we recursively subdivide the cell until the 
weight is acceptably small. Finally, we use the usual rectangle rule and sum over all 
cells to compute the integral. 

It is not really necessary for the map vr : X — )■ to be a multishcetcd cover or even 
to be defined everywhere. In particular, a dominant rational map would be perfectly 
fine. To summarize, our integration algorithm 

• does not require a cell decomposition of the Calabi-Yau manifold, 

• converges as O(^) with the number of points, just like the standard rectangle 
rule, and 

• makes it easy to adapt each integration step to reduce numerical errors. 
5.3 Integrating Over Projective Space 

Thus far, we reduced the integration to one over projective space P" with some conve- 
nient auxiliary metric. We now describe a way of decomposing P" into equal-sized cells 
that is suitable for adaptive subdivision of the cells. First, let us decompose P" into 
n + 1 polydiscs 

n 

P"=UA, Di = [[zo:---:Zi-i:^:Zi+i:---:Zn]\\zj\<l}c,D'\ (5.4) 

*=0 i-th entry 

Note that, by dividing with the homogeneous coordinate of largest magnitude, the 
homogeneous coordinates of any point can be rescaled such that 

• one homogeneous coordinate equals unity, say, Zi = 1, and 

• all other homogeneous coordinates have equal or smaller magnitude \zj\ < 1, j ^ i. 

Hence, a generic point of P" is contained in precisely one of the polydiscs -Dj. The 
polydiscs overlap in the measure-zero set where two or more homogeneous coordinates 
attain the maximum magnitude. 
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coarse 20 points) medium (w 80 points) fine 300 points) 




Figure 3: The adaptive mesh generation for four different integrations over the 
disk D = {\z\ <V\. In each case, the points are chosen such that their 
weight is approximately constant. The point distributions in the three 
rows are further explained in items |7| |^ andj^ on page 
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It remains to decompose each polydisc Di ~ D". For this purpose, we use that the 
polydisc is the Cartesian product of n individual disks 



D = \ re'^ 



r < 1, < (/? < 27r^ C C. (5.5) 



Any cell decomposition of Z) C C induces one of the polydisc. For our implementation, 
we chose to decompose D into annuli of width 5r, and each annulus into segments of 
angles 6(p ~ 5r/(27rr). The annulus segments are approximately quadratic and of area 
2Trr6ip ■ 5r in the Euclidean measure. 

To summarize, we decompose the polydiscs Di into hypercubes of constant volume 
in the Euclidean metric, which we use as the auxiliary measure. If the weight of the cell 
is too large, we subdivide it by splitting the annulus segment in one of the n disc^that 
constitute Di . Summing over each of the n + 1 polydiscs according to the rectangle rule 
then computes the integral over P". To illustrate the adaptive integration algorithm. 



we plot the point distribution for three different volume forms in Figure 3 The first, 
second, and third column shows successive refinements with growing number of points. 
Each row corresponds to a different volume form that the point distribution is adapted 
to: 

1. In the top row, points are distributed regularly on the disc. This distribution, with 
constant weight attached to each point, approximates integration with respect to 
the Euclidean measure. 

2. In the second row, we try to approximate integration with respect to the Fubini- 
Study measure. Here, the points are chosen adaptively and are denser towards the 
center where the Fubini-Study volume form is denser. 

3. In the third row, we illustrate integration over the (2, 2)-hypersurface 

- ^wow.zf - '-^wfzi - m^wUzr + ^-^wlzl = (5.6) 

in Pj^,^.^^] X Pj^^jj.^^] with respect to the Calabi-Yau volume form. Using the pro- 
jection vr : X — 7- Pu(,-2-^] we rewrite this integration as the integration over a single 



^. The disc in Figure 3 shows the first patch [zq : -Zi] = : 1], < 1- The 



points are adaptively chosen to have approximately constant weight. 



6 Diagnosing Stability 

One of the goals of this paper to use the generalized Donaldson algorithm to probe 
Kahler cone substructure in higher- dimensional Kahler cones. That is, we would like to 
efficiently be able to scan a Kahler cone for regions where a given bundle is stable. To 

"^We chose one of the n discs randomly, but this could clearly be improved by separating the cell in 
the direction of the biggest gradient for the weights. 
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this end, we present in this section a new and efficient numerical measure of the stabihty 
properties of V at a fixed polarization. 

The central idea is as follows. If one just wants to numerically determine whether 
a bundle is stable or not, it is not necessary to explicitly compute the Hermitian Yang- 
Mills connection. Instead, one can simply use the first step of the generalized Donaldson 
algorithm, namely the convergence properties of the T-operator. As we saw in |Subse"c^ 
tion 3.2 , according to Wang's theorem, for a fixed embedding defined by H^{X, V^£^"), 



the T-operator can be moved to a balanced place if and only if it is Gieseker stable as 



defined in eq. (3.20). Thus, for an embedding defined by even a relatively small number 



of sections (that is, a small degree of twisting kn), it should be straightforward to check 



whether or not the iterations of the T-operator in (3.19) converge to a fixed point. If 
such a fixed point exists, then we know that the bundle is Giesker stable. In general, 
this computation is easier and faster than a complete computation of the connection 
and its integrated error measure. 
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Figure 4- The T-operator convergence measured for the same bundle and polar- 
izations as in 



Figure 2 



As part of our program of probing the stability of general bundles V in higher di- 
mensional Kahler cones, it is intriguing to see that one can gain important information 
regarding Gieseker stability from the T-operator and only a minimal number of twists 
V (8) C^" . However, one must be careful. While it is certainly true that all slope-stable 
bundles are Gieseker stable, properly Gieseker stable bundles need only be slope semi- 
stable. We can only conclude that if the T-operator fails to converge, then V is not 
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slope-stable. However, if it does converge, we only know that V is semi-stable, but not 



necessarily a solution to (2.1 ). In order to be certain that V is properly slope-stable, one 



would then need to augment the analysis of the T-iteration with a full computation of 



the r error measure in (4.3), as discussed in the previous sections. In the next section. 



we will investigate this subtle semi-stable behavior in detail. For now, however, we only 
explore how much can be gained from a simple check of T-operator convergence. 

With this idea in mind, we develop a new error measure based on the convergence 
of the T-operator. Considering the matrix H^j^ in (3.22), we would like to know how 



much the matrix changes as one moves from the m-th to the (m -|- l)-th iteration of 



the T-operator in (3.19). To answer this question, consider the eigenvalues of the 



matrix. At each step of the iteration, the largest eigenvalues are the relevant features; 
small eigenvalues only give small corrections to the connection on the bundle. Let 



(m) 



(6.1) 



be the largest eigenvalue of the H^^ after m iterations of the T-operator. Except for the 
overall numerical scale, one expects that the details of the initial values are washed out 
by the iteration. However, the overall scale is preserved by the iteration and, moreover, 
does not enter the connection. Therefore, a good quantity to measure the convergence 
is 



(m) 



itut2) 



(m — 1) 



1. 



By Theorem 4 and the above considerations, we know that 



lim r^^^ Jm) 







(6.2) 



(6.3) 



if V is Gieseker stable as in eq. (3.20). For the purposes of this paper, we are interested 
in slope-stability and solutions to the Hermitian Yang-Mills equations. Since slope- 
stability implies Gieseker stability, it follows that '''^^^^^inT-) should have the following 
behavior in the presence of Kahler cone substructure: 



lim r™^^ I 



m) 



= if (tl,t2) G ^stable 
> if (tl,t2) G ^unstable 



(6.4) 



We will put this to the test in the example introduced in Subsection 4.2 For the 



SU{2) bundle V in (4.5), one would expect 



lim r^'^^ ^ 



(m) 



= 
> 



if I 



> 



(6.5) 



The results are shown in Figure 4 for the same six polarizations chosen in Figure 2 , The 



connection T-operator was computed at 30 iterations and the comparison of (6.2) made 
between the 29th and 30th iterations. As expected. 



(30) is approximately zero in 
the slope-stable region of Kahler moduli space, but non-zero in the unstable region. 
On the boundary between /Cgtabie and /Cunstabiei however, we find that the r^^^{?,Q) 



is also close to zero, despite the fact that V in (4.5) is only slope semi-stable for this 
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0.8 




ti:t2 

Figure 5: Convergence of the T-iteration on the SU (2) bundle in ( |4.5[ ) with sub- 
structure for different Kdhler moduli. The radial direction gcd(ti, ^2) is 
always chosen to be as large as possible subject to the constraint that 
there are dim (^0{ti,t2)) < 200 sections. 



polarization. This is to be expected however, since it can be shown via direct computa- 
tion (and (3.20)) that V is still properly Gieseker stable for this line in Kahler moduli 
space and Gieseker stability implies only slope semi-stability. It follows that to accu- 
rately determine the behavior on this boundary, one would need to also investigate the 
r error measure of the previous section, which can distinguish between slope-stable and 
semi-stable behavior. We return to the boundary behavior in the next section. 

For now, it should be noted that one need not have checked the T-operator for all 
500 sections in Figure 4 In fact, the convergence of the T-operator is already evident 
at a much smaller projective embedding. While we must have enough sections to make 
sure that (3.15) is a proper embedding, we can see the stability properties of V from 
the first embedding which makes H'^{X,V C^") non- vanishing. In Figures [s] and [g] 
we present the same results from a purely "angular" point of view; that is, computing 
the T-operator at only a single point along each of the rays plotted in Figure 4 These 
points (that is, the twistings C^^) were chosen so that there are < 200 sections at each 
point. Note that the oscillation in the height of the points in Figure 5 in the unstable 
region is not significant since we are comparing data from different polarizations, which 
leads to different normalizations. The only meaningful comparison is between zero and 



non-zero values of f^^^^-^ 



(m). We also present a 3-dimensional plot of the same Kahler 
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cone substructure in [Figure 6[ For all calculations, the connection was computed with 
Ng = 74, 892 points (adaptive) and 30 iterations of the T-operator. 

Finally, in Figure 7| we consider the rate at which the T-operator approaches its 
fixed point as the number of sections defining the embedding is increased, see eq. (3.15). 
We observe that convergence of the T-iteration is generally slower for more sections. In 
particular, in Figure 7 we present the convergence (and divergence) of the T-iteration 
for polarizations 0{t,t), which is in the unstable region /Cunstabie of the Kahler moduli 
space. 



7 On The Line Of Semi-Stability 



In the previous sections, we demonstrated that the generalized Donaldson algorithm 
is capable of broadly probing Kahler cone substructure. In this section, we take a 
detailed look at the boundary between stable/unstable regions, the so-called "stability 

wall" tsgso]. 



Let us revisit the SU{2) bundle V on the K3 surface defined by (4.5). As discussed 
Subsection 4.2 V is destabilized in part of the Kahler cone by a rank 1 subsheaf 
J" C V defined in ([46|). In the region /Cstabie in ([4?f|), < and in /Cunstabie, 
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Figure 7: 



Convergence (for the polarization 0{2t,t)) and divergence (for 0{t,t)) 
of the T-iteration on the bundle V in (4.5 1. Substructure is shown for 



different Kdhler moduli along a stable and an unstable ray, respectively. 
The orange and gray lines are associated with two different values of 
the bundle moduli for the unstable polarization. 



fj,{J-) > 0. What happens on the boundary hne between these two regions, where 
n{J-) = fJ-C^) = 0? By definition, on a line with /u(J^) = f^{V) the bundle V is semi- 
stable. Hence, for generic values of the bundle moduli its connection will not solve the 
Hermitian Yang-Mills equations. However, semi-stable bundles are distinguished from 
unstable bundles in that they can provide supersymmetric solutions for special loci in 



(2.9 


) and ( 


2.10 


) in 



the only way that V can satisfy the Hermitian Yang-Mills equations when n{J-) = A'(V) 
is for it to be poly-stable rather than strictly semi-stable. That is, if the connection 



on "V is decomposable, it is possible that a connection may exist which satisfies (2.1). 
Mathematically, this property of slope semi-stable bundles is characterized by the notion 
of S-equivalence classes and the Harder-Narasimhan filtration [42i|45j, which states that 
every semi-stable bundle has a unique poly-stable representative in its moduli space. 
For the bundle in (4.5), we find that the poly-stable representative arises when V 



decomposes on the semi-stable wall as the direct sum 

V^^eO(-2,l). (7.1) 
One can quantify this decomposable locus in the bundle moduli space as follows. The 
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Figure 8: The t error measure associated with the SU(2) bundle in (4.5) at the 
stability wall in Kdhler moduli space. The bundle modulus ip defined 
in (7.2) determines whether V is strictly semi-stable (when if ^ 0) or 
poly-stable (when (p — 0). In the latter case, the reducible connection 
satisfies the Hermitian Yang-Mills equations (2.2 1. 



bundle moduli space of V is described by the parameters in the map / in (4.5). We 



can describe the decomposable locus in (7.1) in terms of this map by parametrizing one 



particular vector bundle modulus as 99 in the monad map 



(9 - i)yo + (-6 - 6i)yi + (-8 - 8i)y2 
(-9 - 9i)yo + (-2 + 9i)yi + (4 - 4i)y2 

-10+7j)x;^+(6-9i)xgxi+(3-8j)x§x2+(-9-4i)a;oa;f+(8+4i)a;f)^ 



which determines V, 



(7.2) 







0(-2,-l) 



O(-2,0)®2 + 0(2,-1) 







(7.3) 



When (f = 0, the bundle splits as in ( |7.1| ) and the resulting direct sum is poly-stable and 
a solution to the Hermitian Yang-Mills equations. When 7^ 0, V is only semi-stable, 



one cannot solve (2.2), and supersymmetry is broken. 



Let us now revisit the numerical results of the past few sections in light of this 
structure. Can the error measures ri^Ay) and r™'"^(m), introduced in Subsection 4.1 



and Section 6 respectively, accurately reveal this (/9-dependent structure? To answer this 
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Figure 9: The convergence/divergence of T-iteration associated with the SU{2) 



bundle in (4.5) at the stability wall in Kdhler moduli space. The bundle 



modulus If defined in (7.2) determines whether V is strictly semi-stable 
(when ip ^ 0) or poly-stable (when = 0). In the latter case, the 



reducible connection satisfies the Hermitian Yang-Mills equations (2.2) 



question, we have run the algorithm for the same bundle again, this time keeping the 
Kahler moduli fixed directly on the line (ti, ~ (4, 3), but allowing the bundle moduli 
to vary through the variable (p. The metric was computed for this polarization with 
Ng = 202, 800 points (adaptive), 30 iterations of the metric T-operator and with a Kahler 
form determined by 0(4, 3)®^. Meanwhile, the connection utilized Nq = 74, 892 points 
(adaptive) and 100 iterations of the connection T-operator. The results for the r(ii,t2) 
error measure of equation (4.3) are shown in Figure 8 Significantly, the algorithm 



clearly distinguishes between the semi-stable {(p ^ 0) and poly-stable (93 = 0) behavior 
of V. Moreover, one can repeat the same analysis with the T-operator convergence 



measure, r"^"^(12, 9)(m), of (6.2). For this analysis, we hold ourselves fixed at one point 
in Kahler moduli space (ti = 12, t2 = 9) and allowing the number of T-iterations m to 
increase. Once again, we find that this error measure definitively distinguishes between 
the supersymmetric and non-supersymmetric configurations. The convergence of the 



T-iteration is shown in Figure 9 



The graph [Figure 9 shows a high degree of sensitivity to the bundle moduli. The 
fact that the T-iterations can distinguish the </? = and ip ^ cases is significant, 
since it means that for a general value of the bundle moduli V is not Gieseker stable 



30 



for the boundary wall polarization. This can be verified directly using the destabilizing 
subsheaf J- in (4.6) and the definitions of Gieseker stability, (3.20), in Section 2 In fact, 
for £ = 0(4, 3) we have 



while 



Therefore, 



P£(-F)(n) =45n"-3, 
P£(V)(n) = J(90n2 -8) 



Pz:(V)(n) -<p£(-F)(n) 



(7.4) 
(7.5) 
(7.6) 



and, by (3.21), V is not Gieseker stable. The results of Figure 9 are in complete agree^ 



ment with this; showing the = reducible connection T-operator converging to its 
fixed point while the Gieseker unstable configurations with 99 7^ fail to converge. It fol- 
lows that the T-operator is sensitive enough to the bundle moduli that it can distinguish 
between these important cases. 



8 A Calabi-Yau Threefold Example 

In the previous sections, we investigated Kahler cone substructure with a number of 
new tools and from a variety of perspectives. Here, we present a final example involving 
geometry directly relevant to realistic M = \ heterotic compactifications; namely, a 
Calabi-Yau threefold with a higher-rank vector bundle defined over it. 

As a base manifold, we choose a Calabi-Yau threefold determined by a generic degree 
(4, 2)-hypersurface in p3 X For this threefold, /i^'^ = 2 and all line bundles can be 
written in the form 0{ai, 02), where 0{1, 0) arises from the hyperplane descending from 

and 0(0, 1) descends from the hyperplane of P^. Over this threefold, define the SU{3) 
monad bundle 

0^ V^O(-2,l)©O(3,-l)eO(2,0)eO(2,l)®2 ^0(4,l)eO(3,l) ^0 (8.1) 
and its associated dual bundle 

0(-4, -1) © 0(-3, -1) 

0{-2, -1)®2 © 0{-2, 0) © 0(-3, 1) © 0(2, -1) ^ ^ 0. (8.2) 

This example was chosen because it once again the stability depends on the precise value 



of the Kahler moduli. It can be verified, using the definitions of Section 2 , that the bundle 



V in eq. (8.1) has two potentially destabilizing sub-sheaves. These are Ti,T2 C V given 
by 

^ -Fi ^ 0(3, -1) © 0(2, 0) © 0(2, 1)®2 Jl^ 0(4, 1) © 0(3, 1) ^ 0, 

(8.3) 

J-2 ^ 0{-2, 1) © 0{2, 0) © 0{2, 1)®2 ^ 0(4, 1) © 0(3, 1) 0. 
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Figure 10: 



Convergence of the T-iteration on the SU(3) bundle defined by (8.2 1 
for different Kdhler moduli. The radial direction gcd(ti,f2) is always 
chosen to be as large as possible subject to the constraint that there are 
dimi7'^(0(ii, ^2)) < 500 sections. The results agree with the Kdhler 
cone substructure in (8.4). 



Both are of rank 2 with ci(J-'i) = (2,-1) and ci(J-2) = (—3,1). Due to these two 
destabilizing sub-sheaves, the Kahler cone divides into three regions 

/Cstable = {0(il,i2) I I < I < |}, 
/Cunstable = {0(tl,t2) || < |} U {0(il,t2) || > |}- 

Moreover, a direct calculation using ( |3.20| ) shows that V is still Gieseker stable at each 
of the stability walls defined by t2/ti = 4/3 and t2/ti = 4/5. Thus, while the bundle 
is strictly semi-stable along these rays (and, hence, does not solve the Hermitian Yang- 



Mills equations eq. (2.2)), one would still expect the T-iteration to converge for points 
along these boundaries. 

For ease of embedding, the T-operator is computed for V^. However, since V and 
its dual must be slope-stable/unstable in the same regions of Kahler moduli space, this 
does not affect our results. Since we have demonstrated in the previous sections that the 
T-iteration is a fast and efficient check of stability, we will consider here only the error 
measure '"(^"tj)^™)' interesting to observe that, even in this more complex example, 
the data shown in [Figure 10| clearly reproduces the Kahler cone substructure described 
in (8.4). Moreover, at the points on the two stability walls we find, as expected, that 
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r™^''^2)('^^) ~ ^' since V is still Gieseker stable for these lines of slope semi-stability. 
The connection integration was performed with Nq = 119, 164 points (adaptive) and 30 
iterations of the T-operator. The plot shown in Figure 10 are analogous to the results 
in 



Figure 5 for the K3 surface. 



9 Conclusions and Future Work 

In this paper, we extended the generalized Donaldson algorithm to manifolds with 
higher-dimensional Kahler cones and presented a method for approximating the con- 
nection on slope-stable holomorphic vector bundles in this context. We also introduced 
a numerical criterion for determining the existence of supersymmetric heterotic vacua, 
without having to compute the connection. These techniques clearly can be used to 
search for new classes of smooth M = 1 supersymmetric vacua in heterotic string the- 
ory. However, the explicit knowledge of the gauge connection that they provide allows 
us to go far beyond a simple categorization of vacua. 

It is a long-standing goal of string theory to produce low-energy theories whose effec- 
tive actions reproduce the symmetries, particle spectrum, and properties of elementary 
particle physics. Within the context of smooth heterotic string or M-theory compactifi- 
cations [21 [71 160] , there has been considerable progress in recent years [611 [62l [631 EH IMl 
[nni[S3]- Despite these successes, certain observable quantities of particle physics, such as 
the gauge and Yukawa couplings, are difficult to compute directly. Normalized Yukawa 
couplings, for example, depend on both the coefficients of cubic terms in the superpo- 
tential and the explicit form of the Kahler potential. In turn, these quantities depend on 
the detailed structure of the underlying geometry - that is, the metric and gauge con- 
nection on the Calabi-Yau threefold and holomorphic vector bundle respectively. Hence, 
the Yukawa couplings in the four-dimensional effective theory are not known, except in 
very special cases. One can rarely do better than the qualitative statement that such 
coupling coefficents either vanish or are of order one. However, the methods introduced 
in |32[l29j and extended in this paper allow one, for the first time, to explicitly compute 
both the metric and the gauge connection and, hence, the Yukawa coefficients. In future 
work [HT], we will use these techniques to calculate these physical couplings. 
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